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ABSTRACT 

Context. The star formation rates in starburst galaxies are orders of magnitude higher than in local star-forming regions, and the 
origin of this difference is not well understood. 

Aims. We use sub-mm spectral line maps to characterize the physical conditions of the molecular gas in the luminous Galactic star- 
forming region W49A and compare them with the conditions in starburst galaxies. 

Methods. We probe the temperature and density structure of W49A using H2CO and HCN line ratios over a 2' x 2' (6.6 x 6.6 pc) 
field with an angular resolution of 15" (~0.8 pc) provided by the JCMT Spectral Legacy Survey. We analyze the rotation diagrams of 
lines with multiple transitions with corrections for optical depth and beam dilution, and estimate excitation temperatures and column 
densities. 

Results. Comparing the observed line intensity ratios with non-UTE radiative transfer models, our results reveal an extended region 
(about r X r, equivalent to ~3 x 3 pc at the distance of W49A) of warm (>100 K) and dense (>10^ cm"') molecular gas, with a 
mass of 2 X lO'' - 2 x 10^ Mq (by applying abundances derived for other regions of massive star-formation). These temperatures and 
densities in W49A are comparable to those found in clouds near the center of the Milky Way and in starburst galaxies. 
Conclusions. The highly excited gas is likely to be heated via shocks from the stellar winds of embedded, O-type stars or alternatively 
due to UV irradiation, or possibly a combination of these two processes. Cosmic rays. X-ray irradiation and gas-grain coUisional 
heating are less likely to be the source of the heating in the case of W49A. 

Key words, stars: formation - ISM: molecules - ISM: individual objects:W49A 

1 . Introduction HCO^ J=l-0 spectra measured toward the ring-like configura- 
tion of UC Hii regions in terms of an infall profile. Serabyn et 
The formation of stars proceeds on scales ranging from 1000 al. (1993) interpret the two-components of the double-peak line 
AU-sized isolated low-mass cores to -10 pc-sized GMCs profile seen in CS (J=10-9, 7-6, 5-4, 3-2) and C^'^S (J=10-9, 
with embedded massive clusters. The observed levels of star- 7_6, 5-4) transitions as coming from different clouds, and sug- 
formation activity also vary gready: from -10"^ Mg/yr in dwarf gest that massive star formation in W49A is triggered by a large- 
galaxies like I Zw 18 (Aloisi et al. 1999) to ~3 Mo/yr in the scale cloud merger. Williams et al. (2004) argue that one cannot 
Milky Way to -10^ Mo/yr in starburst galaxies (e.g. Solomon distinguish between a global collapse model of a large cloud and 
. & Vanden Bout 2005). The origin of this large spread in star a multiple-cloud model, based on HCO+ (J=3-2 and 1-0) and 
; formation rate may lie in different physical conditions with the c'**0 (J=2-l) data. Roberts et al. (201 1) observe single-peaked 
! gas forming the stars, the nature and role of feedback effects or HC'^0+ (J=4-3) line profiles, which supports the large-scale in- 
H • other mechanisms. As one of the the most active star-forming fall hypothesis. In '^^CO and C"^0 J=2-l maps of the W49A re- 
? : regions in the Galactic disk, W49A provides a relatively nearby gion Peng et al. (2010) find evidence of expanding shells with the 
laboratory in which to study intense star formation activity and dense clumps in the region lying peripherally along the shells, 
may provide insights into explaining the large observed range in These expanding shells could provide a natural explanation of 
star-formation activity in galaxies. the observed molecular line profiles, and triggered massive star 

W49A is a luminous (>10^ Lq) and massive (-10^ Mq) star- formation in W49A. 
forming region (Sievers et al. 1991) at a distance of 11.4 kpc 
(Gwinn et al. 1992). W49A contains several signatures of high 
star-formation activity, in particular, strong H2O maser emission 

(Genzel et al. 1978). Line ratios of HNC, HCN and HCO+ and xhis paper reports new kinetic temperature and density es- 

theirisotopologues suggest similar conditions to starburst galax- timates for a 2' x 2' region of W49A based on single dish, 

ies (Roberts et al. 2011). Several ideas have been proposed to jcMT images with 15" resolution. These indicate the presence 

explain the high star-formation activity of W49A. of extended warm and dense gas towards the central region of 

Welch et al. (1987) propose a large-scale inside-out grav- the source. We use these results to compare W49A to starburst 

itational collapse by interpreting two velocity components of galaxies where similar results have been reported. 
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2. Observations and Data reduction 

This paper presents results from data obtained as part of the 
JCMT Spectral Legacy Survey (SLS, Plume et al. 2007) cov- 
ering between 330 and 373 GHz which has been carried out at 
the James Clerk Maxwell Telescope (JCMT)' on Mauna Kea, 
Hawai'i. The data were taken using the HARP-B receiver and 
the ACSIS correlator (Buckle et al. 2009). HARP-B consists of 
16 pixels or receptors separated by 30 arcsec with a foot print of 
2 arcmin. 

The observations were carried out using jiggle - position 
switch mode. Each map covers a 2' x 2' field sampled every 
7.5", half of the JCMT beam width at the relevant frequencies 
(15" at 345 GHz). The spectral resolution at these frequencies is 
~0.8 km s^'and the beam efficiency is 0.63 (Buckle et al. 2009). 
The spectra were taken using an off-position 14' to the northeast 
of the source. 

Data reduction was done manually using the tools in the 
Starlink software package. The raw time series files are pro- 
cessed in 1 GHz wide chunks. These files were inspected for 
bad receptors, baselines and frequency spikes. Bad receptors, 
bad baselines and spikes were masked in the time series files be- 
fore they have been converted to three-dimensional data cubes 
(RAxDecxfrequency). These 3D cubes then were baseline cor- 
rected by fitting and removing a linear or second-order polyno- 
mial baseline. The baseline corrected cubes are then mosaicked 
to a larger cube covering an area of 2' x 2' and the frequency 
range of interest. Single sideband system temperatures for these 
observations are in the range between 255 K (at 360 GHz) and 
700 K (at 373 GHz). The analysis here primarily focuses on a 
number of spectral lines in the 360-373 GHz range plus a few 
selected lines from the 330-360 GHz part of the survey. The 
whole, systematic data analysis of the full SLS frequency range 
is the focus of a future paper. 

We also use ancillary observations of the HCN J - 3-2 
transition made using JCMT raster maps with half-beamwidth 
spacing using the Receiver A (211-276 GHz; beamwidth 20"; 
main beam efficiency: 0.69). The HCN 3-2 map covers approx- 
imately the same field that of the SLS HARP maps. The HCN 
3-2 data were reduced and calibrated with standard Starlink 
procedures. Linear baselines were removed and the data were 
Hanning smoothed. The map is a grid of 15 pixel x 15 pixels, 
where each pixel covers 8"x 8". 



3. Results 

3. 1 . Line selection 

Figure 1 shows the 360-373 GHz spectrum towards the central 
position (RA(J2000)=19:10:13.4; Dec(J2000)=09:06:14). More 
than 80 lines have been identified in this spectrum. 

Our initial selection of lines for the analysis here is based on 
molecules that show multiple transitions in the high-frequency 
range (360-373 GHz) of the SLS, which include six H2CO, five 
CH3OH and nine SO2 transitions (Table 1). These are used to 
trace the excitation conditions. One more H2CO transition from 
the lower frequency part of the survey (H2CO 5 15 - 4 14) was 
added to be used as a tracer of kinetic temperatures, together 
with the 533 - 432 transition from the high-frequency part of the 



' The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 



Table 1. Summary of the molecular lines used in this paper 
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" noise level measured over the central 1' x 1.5', for 1 km s ' chan- 
nels. 



SLS. The HCN 4-3 transition was selected from the lower fre- 
quency range, to trace the density structure, together with the 
HCN 3-2 transition. 

In this paper, we restrict our analysis to our line selection 
summarized above and in Table 1. A subsequent paper (Nagy 
et al, in prep.) will give more details on the chemical inven- 
tory including the full SLS frequency range (330-373 GHz) 
and on kinematics traced by various species with significant 
spatial extension. Figure 2 shows the main regions around the 
center of W49A overplotted on the integrated intensity maps 
in the HCN 4-3, SO2 152,i4 - 14i,i3, CH3OH 72,5 - 61,5 and 
H2CO 5 15 - 4i4 lines. The integrated intensities of the HCN 
and H2CO lines show a similar distribution, indicating that 
their emission originates in the same volume. This is also sup- 
ported by the line profiles of these two species (Nagy et al, 
in prep). In the following sections, we will refer to these re- 
gions with the coordinates indicated in the figure: source center 
(RA(J2000) = 19:10:13.4; Dec(J2000) = 09:06:14), Eastern tail 
(RA(J2000) = 19:10:16.6; Dec(J2000) = 09:05:48), Northern 
clump (RA(J2000) = 19:10:13.6; Dec(J2000) = 09:06:48) 
and South-West clump (RA(J2000)=19:10:10.6; Dec(J2000) = 
09:05:18). 



3.2. Excitation conditions 

A well known method for the estimation of excitation tempera- 
tures and column densities is the rotational diagram method (e.g. 
Turner 1991). It can be applied to molecules with multiple ob- 
served transitions using three assumptions: 1) the lines are op- 
tically thin, 2) the level populations can be characterized by a 
single excitation temperature ('rotational temperature', Tiot) and 
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Fig.l. The spectrum towards the source center (RA(J2000)=19:10:13.4; Dec(J2000)=09:06:14) in the 360-373 GHz frequency 
range of the SLS. The main detected Unes are labeled. The gap in the spectrum around 368.5 GHz is a region of high atmospheric 
extinction and was not observed. 



3) the emission is homogeneous and fills the telescope beam. 
Then the measured integrated main-beam temperatures of lines 
(/ Tmb^^V K km s ') is related to the column densities of the 
molecules in the upper level (A^u) by: 

Au Mot / EA 1.67 -lO'^ r 

— = exp U — — TMBdV, (1) 

with gu the statistical weight of level u, A^tot the total column 
density in cm"-, Q{Ti-ot) the partition function for Tjot, the 
upper level energy in K, v the frequency in GHz, yU the permanent 
dipole moment in Debye and S the line strength value. A linear 
fit to ln(Au/gu) - Ea gives Ti-ot as the inverse of the slope, and Atot, 
the column density can be derived. The rotational temperature 
would be expected to be equal to the kinetic temperatures if all 
levels were thermalized. 

In the case of non-LTE excitation, a different excitation tem- 
perature may characterize the population of each level relative to 
that of the ground state or relative to that of any other level. In 
this case, the three assumptions mentioned for the rotation dia- 
gram method do not apply. Equation 1 can be modified to include 
the effect of optical depth t through the factor Cr=T/(l-e"^) and 
beam dilution / (= Qi/0„) (Goldsmith & Langer 1999), with Q, 
the size of the emission region and 0„ the size of the telescope 
beam. 

In ( ^) = In ( - + In (/) - In (C,). (2) 

According to equation 2, for a given upper level, N°^^ can be 
evaluated from a set of Atot.i, T^ex, / and Cj- Since Cj is a func- 
tion of Atot,i and Tex, the independent parameters are therefore 
Ntot,\ , Tex and /. Solving the equation for a 'source size' between 
1" and 15" (which is equivalent with a beam dilution factor be- 
tween (1/15)^ and uniform beam filling), minimization gives 
Atot.i and Tgx. Figure 3 shows the rotation diagrams for CH3OH, 
H2CO & SO2 toward the central position with (green symbols) 
and without (red symbols) taking into account the optical depth 
and beam dilution. Table 2 includes the parameters of the rota- 
tion diagram fit toward the central position. 

The CH3OH lines are found to be optically thin (r = 0.004 - 
0.02), so the rotation diagram fit is a good approximation to the 
column densities and excitation temperatures. The derived rota- 
tion temperature is 83+7 K. The H2CO and SO2 lines are mostly 
optically thick with optical depths in the range of 0.2-2.35 for 



H2CO and 0.1-5.4 for SO2. The best fit excitation temperatures 
are 138 K for H2CO and 1 15 K for SO2. The best fit total column 
densities are 4.3 x 10"' cm"- forH2CO and 1.4 x 10'** cm"- for 
SO2. The;^^- minimization for H2CO and SO2 results in a best fit 
'source size' of ~2-3", which is equivalent with 0.11-0.17 pc 
for W49A. This is consistent with the sizes of hot cores and UC 
Hii regions resolved by sub-arcsecond resolution measurements, 
such as Wilner et al. (2001), De Free et al. (1997, 2004). 

Figure 4 shows the rotation diagrams for H2CO toward the 
three main regions outside the center The rotation diagram 
with corrections for optical depth and beam dilution (Figure 
4. and Table 3.) gives similar excitation conditions toward the 
three main sub-regions (excitation temperatures; 67 K for the 
Northern clump, 56 K for the Eastern tail and 6 1 K for the South- 
west clump). These excitation temperatures are about half of 
the excitation temperature towards the central position (138 K). 
Since the lines are optically thin, the rotation temperatures are 
good estimates for the excitation temperatures. The best fit col- 
umn densities are also similar between the off-center positions, 
about 1.7-2.2x10'"' cm"-', almost three orders of magnitude less 
than the best fit column density towards the center The data at 
these positions are consistent with optically thin emission and 
uniform beam filling. 

3.3. Kinetic temperature estimates 

To estimate the local physical conditions in W49A, we use ratios 
of lines from the same molecular species. The kinetic tempera- 
ture can be probed by the line ratios of H2CO lines from the same 
y-state but different /T-states. In this case, we use the 0-H2CO 
5 15 - 4i4 (35 1 .769 GHz) and 0-H2CO 533 - 432 (364.275 GHz) 
transitions, since they have been detected in a large part of the 
SLS field and have been found to be excellent tracers of kinetic 
temperature (Figure A.l). Our calculations use the non-LTE ra- 
diative transfer program Radex (Van der Tak et al. 2007), molec- 
ular data file from the LAMDA database (Schoier et al. 2005) 
and H2CO collision data from Green (1991), which have been 
scaled by 1.37 to make a first order approximation for collisions 
with H2. Uncertainties introduced by using these rates may be up 
to 50% for the column densities, but less affect the densities and 
temperatures, since those are derived from line intensity ratios. 
More recent calculations by Troscompt et al. (2009) are more 
accurate, however, include a more limited range of temperatures 
and transitions, therefore, were not used in our calculations. 
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Fig. 3. The results of the rotation diagram analysis for H2CO to- 
ward the center Top panels: The results of the rotation diagram 
analysis before (red symbols) and after (green symbols) correc- 
tions for optical depth and beam dilution. The overplotted blue 
line corresponds to a linear fit to the rotational diagram without 
corrections. Bottom panels: The corresponding best fit optical 
depths from the minimization. 
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Fig. 2. The integrated intensity distribution of HCN 4-3, SO2 
152,14- 14i, 13, CH3OH 72,5-61,5 and H2CO 5i5 -4u- The central 
position is shown with black plus signs, the off-center regions 
(Northern clump. Eastern tail and South-west clump) are also 
shown with white plus signs. 
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Fig. 4. The results of the rotation diagram analysis for H2CO at 
off-center positions. Top panels: The results of the rotation di- 
agram analysis before (red symbols) and after (green symbols) 
corrections for optical depth and beam dilution. The overplot- 
ted blue line corresponds to a linear fit to the rotational diagram 
without corrections. Bottom panels: The corresponding best fit 
optical depths from the x~ minimization. 



Figure A.l shows the calculated line ratios as a function of 
kinetic temperature and H2 density. In the central 40" x 40" re- 
gion, FWHM=12 km s"' was used as an initial parameter for 
Radex, while outside of this region we used FWHM=6 km s 
based on the observed values. However, there is no strong depen- 
dence on the FWHM of the lines. We apply a molecular column 
density of A^(H2CO)=5 x lO'"* cm"^ as an 'average' value for all 
the positions over the SLS field. This is slightly below the value 
that can reproduce the optical depths derived from rotation di- 
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agram method for the central position, assuming uniform beam 
filling. 

The results are shown in Figure 5 with the main regions high- 
lighted. We derive a kinetic temperature of 300 K toward the 
center, 98 K toward the Northern clump, 85 K toward the South- 
west clump and 68 K toward the Eastern tail, with an uncertainty 
of ~30%. In fact, there is a large, about 1' x 1' (~ 3 x 3 pc) region 
around the center, with kinetic temperatures >100 K. 

3.4. Volume- and column density estimates 

To estimate H2 volume densities of the high temperature gas that 
the H2CO line ratios trace, we use line ratios of the HCN 3-2 
(265.886 GHz) and HCN 4-3 (354.505 GHz) transitions. Our 
calculations are based on collisional rates from Dumouchel et 
al. (2010), which have been scaled by a factor 1.37 to represent 
collisions with H2. The HCN 4-3 map has been convolved to the 
resolution of the HCN 3-2 map (20")- Since the line ratios de- 
pend on the kinetic temperature (Figure A. 2), we use the kinetic 
temperatures determined from the H2CO line ratios to convert 
the measured HCN line ratios (Figure 6) to volume densities. 
Therefore, even if HCN was detected towards almost the entire 
SLS field, our calculations only take into account the positions 
where both of the H2CO transitions were detected. The uncer- 
tainty of these density estimates is ~ a factor of two. Figure 7 
shows the estimated volume densities. Note that not all the posi- 
tions with kinetic temperature estimates are included in the den- 
sity map: this is due to the uncertainty of our method below 100 
K, where the HCN 3-2 / 4-3 line ratios are rather temperature 
than density tracers. With a search range in the densities between 
10^ and 10^ cm"^, an extended area towards the center (approxi- 
mately corresponding to the high-excitation area revealed by the 
H2CO line ratios) corresponds to densities as higher than 10^ 
cm"^ or higher, with a maximum of 6.9 x 10*" cm""*. The max- 
imum density is measured toward an offset position compared 
to the position we refer to as 'center'. The densities toward the 
central part show little less than an order of magnitude variation. 

The average H2 volume density corresponding to the central 
1x1 arcminute region is 1.6 x 10* cm"^. The lowest densities, 
toward the edges of the cloud area covered by our estimates, 
is a few times 10"* cm"-'. However, our density estimates are 
more accurate for densities between 10^ and 10^ cm"^^ (Figure 
A. 2). In addition to the uncertainty of the method below ~10^ 
cm"-', the H2CO 5i5 - 4i4 and 533 - 432 transitions are detected 
over a less extended region, compared to the HCN 3-2 and 4-3 
transitions. Estimating the physical conditions outside the area 
covered by our kinetic temperature and volume density maps re- 
quires molecular line tracers with a larger spatial extent. 

Based on the kinetic temperatures and volume densities pre- 
sented above, we estimate the column densities of the two main 
molecular tracers used in this paper, HCN and H2CO. These es- 
timates use H2 volume densities of 1.8 x 10* cm"-' for the cen- 
ter, 1.2 X 10* cm"-' for the Eastern tail, 5.6 x 10^ cm"-' for the 
Northern clump and 1.3 x 10^ cm"^^ for the South-west clump. 
The kinetic temperatures used for each position are listed in sec- 
tion 3.3. To calculate the HCN column density toward the center, 
we use the H'^^CN 4-3 (345.3 GHz) transition, and ^-C/^^C^ll 
(Wilson & Rood 1994). HCN peaks toward the center, with a 
column density of 7.3 x lO'^^ cm"^. The HCN column densi- 
ties of the Eastern tail and Northern clump regions are com- 
parable: A^(HCN) = 7 - 8 X 10'^ cm"^ toward the Eastern 
tail and A^(HCN) = 1.5 - 2 x lO"^ cm"^ toward the Northern 
clump. N(HCN) toward the South-west clump is a factor of 4-5 
above the column density of the Northern clump: A^(HCN) = 



HCN 4-3 / HCN 3-2 




50 -50 

ARA ["] 



Fig. 6. Map of the line ratio of HCN J = 4-3 to HCN / = 3-2. 
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Fig. 7. Volume density estimates from the HCN 3-2 and HCN 
4-3 lines. 



6.2 - 9.6 X 10"* cm"2. NiHzCO) is -3.2x10'^ - 4.5 x lO'"* cm"^ 
toward the center as well as the South-west clump. The Northern 
clump has A^(H2CO) slightly below the center and South-west 
clump,7.0xl0'^-l.lxl0''* cm , while toward the Eastern tail, 
we measure a column density of 4.9 x 10'-' -5.0 x 10'-' cm"^. The 
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Fig. 5. Left panel: Ratio of 533 - 432 to the 5 15 - 4 14 H2CO transition. Right panel: Kinetic temperature map of W49A determined 
from line ratio shown in the left panel. 



two orders of magnitude difference in the H2CO column density 
toward the center compared to the best fit column density de- 
rived from the rotation diagrams can be accounted for the simple 
assumption used in the rotation diagram method, that the level 
populations can be described by a single excitation temperature. 

The mass of the warm and dense gas can be estimated based 
on the H2CO column density, by assuming an H2CO abundance. 
Van der Tak et al. (2000) derive H2CO abundances in the range 
of lO^'^-lO"'" for thirteen regions of massive star formation. For 
an 'average' column density of lO^'* cm"^ and an abundance of 
10"^, the mass of the central 1' x 1' region with warm and dense 
gas is ~2xl0'* Mq, which becomes ~2xl0^ Mq for an H2CO 
abundance of 10"'°. These are comparable to the mass esti- 
mates based on dust continuum measurements. Ward-Thompson 
& Robson (1990) estimate a mass of ~2.4xl0'' Mq correspond- 
ing to the central ~2 pc. Buckley & Ward-Thompson (1996) 
adopt a value of ~ 10^ Mq as a best fit based on dust contin- 
uum as well as on CS data published by Serabyn et al. (1993). 

A second estimate of the mass of dense gas can be obtained 
from the HCN analysis. Adopting an average number density of 
1.6 X lO*" cm""* over the central 1' x 1' region and assuming an 
approximately spherical region would imply a mass of ~2xl0'' 
Mq. This is a factor of 10-100 larger than the H2CO suggests 
indicating that the dense HCN material fills a fraction between 
1% and 10% of the volume of this extended region. 

Further evidence for a small volume filling factor is given 
by the best fit source sizes from the rotation diagrams toward 
the central position. The column densities derived from Radex 
are beam-averaged, while the best fit column densities from 
the rotation diagram correspond to clumps with sizes indicated 
by the best fit source size. Comparing the HCN column densi- 
ties to the H2 volume densities gives another indication of the 
volume filling factor. While the volume densities change be- 
tween 10^ and a few times 10^ cm"^ and show a rather uniform 
distribution, A^(HCN) shows a 2 orders of magnitude variation 



(7xl0'''-7xl0'^ cm"^) and peaks toward the center. The emis- 
sion may originate in warm and dense gas in a number of unre- 
solved clumps. Such a behavior was found before, by Snell et al. 
(1984) and Mundy et al. (1986) for M17, S140 and NGC 2024. 



Table 2. Results of the excitation analysis toward the center of 
W49A 





I. Rotation diagrams 




SO2 CH3OH 


H2CO 


r,-o. [K] 


156+11 83±7 


182+84 




14.56+2.99 6.58±1.38 


0.78±0.17 


Source size["] 


■2 1 +0.0 1 9 ir+2.5 




Tex [K] 


115!^o 82.1+[6[ 


138.2!l|^;g 


Mou* [10'5 cm-2] 


ituu.u ='■■'-3.2 


4^ + 1857 




II. Non-LTE calculation 




300 




A'LHzCO] [cm-2] 


3.2-4.5x10'* 




A'[HCN] [cm--] 


7.3x10'^ 






LSxlO** 





" Total beam-averaged column density calculated using the assump- 
tion that the lines are optically thin, the level populations can be char- 
acterized by a single excitation temperature and the emission fills the 
telescope beam. 

Total source-averaged column density given by the x' minimiza- 
tion, including a coiTection for the optical depth and a beam dilution 
factor. The errors show the 1-cr upper and lower limits of the best fit. 

'' The uncertainty of the kinetic temperatures is ~30%. 

'' The uncertainty of the density estimates is a factor of 2. 
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Table 3. Results of the excitation analysis toward the off-center 
regions of W49A 





I. Rotation diagrams 




Northern 


Eastern tail 


South-west 




clump 




clump 


r„„ [Ki 

lot L^*-J 


67+13 


56+3 


62±11 


Mo," [10'3 cm-2] 


2.26±0.48 


1.82±0.46 


1.68±0.36 


Source size["] 


14.44+°f^ 


14.44!°f^ 


1 4 44+0.56 


T'cx [K] 


67+-' 

"'-13 


55!^^ 


61+" 


A'toi.i [10 cm J 


9 9+44.8 
^•^-0.2 


1 Q+39.1 
^■'-0.3 


J 7+22.3 
^■'-0.2 




II. Non-LTE calculation 




Northern 


Eastern tail 


South-west 




clump 




clump 


Tun [K]'- 


98 


68 


85 


A'[H2CO] [cm--] 


0.7-1.1x10"* 


4.9-5.0x10" 


3.2-4.5x10"* 


^[HCN] [cm-2] 


1.5-2x10"* 


7-8x10" 


6.2-9.6x10"* 


n [cm-^Y 


5.6x10^ 


1.2x10* 


1.3x10' 



" Total beam-averaged column density calculated using the assump- 
tion that the lines are optically thin, the level populations can be char- 
acterized by a single excitation temperature and the emission fills the 
telescope beam. 

* Total source-averaged column density given by the minimiza- 
tion, including a coixection for the optical depth and a beam dilution 
factor. The eiTors show the 1-tr upper and lower limits of the best fit. 

The uncertainty of the kinetic temperatures is ~30%. 

The uncertainty of the density estimates is a factor of 2. 



4. Discussion 

We have detected warm (>100K) and dense (>10^ cm-^) gas ex- 
tended over a 1' X 1' (~3x3 pc) field towards W49A. Based on a 
rotation diagram analysis with corrections for optical depth and 
beam dilution, we derive excitation temperatures of 115 K (for 
SO2), 82 K (for CH3OH), 138 K (for H2CO) towards the center 
From H2CO line ratios with the same J and different /T-states, 
we derive kinetic temperatures of ~300 K toward the center, 68 
K toward the Eastern tail, 98 K toward the Northern clump and 
85 K toward the South-West clump. Our kinetic temperature es- 
timates are consistent with the excitation temperatures derived 
from the rotation diagram method, given that non-LTE condi- 
tions are expected, with excitation temperatures below the ki- 
netic temperatures. Ratios of HCN 3-2 / HCN 4-3 line intensi- 
ties, using our kinetic temperature estimates, result in H2 volume 
densities of 1 .8 x 10'' cm""' for the center, 1 .2 x 10* cm""* for the 
Eastern tail, 5.6 x 10^ cm"^ for the Northern clump and 1 .3 x 10^ 
cm"^ for the South-west clump. 

Our kinetic temperature estimates are significantly higher 
than previously reported values. Based on Planck brightness 
temperatures derived from the CO 7 = 7-6 transition, Jaffe et 
al. (1987) estimate kinetic temperatures of >70 K in the centre, 
falling below 50 K in the outer parts of the cloud. However, these 
observations have a beam size of 32", which suggests higher ki- 
netic temperatures on a smaller scale, and as such, may be con- 
sistent with our results. Roberts et al. (201 1) use HCN/HNC line 
ratios compared to chemical models to estimate kinetic temper- 
atures and derive temperatures of 100 K toward the center, 40 K 
toward the South-west clump, 40 K toward the Eastern tail and 
75 K toward the Northern clump. One possible reason of the dif- 
ference between the kinetic temperatures presented in this paper 
and the values estimated by Roberts et al. (201 1) is that Roberts 



et al. (2011) adopted the HCN colHsion rates for HNC, which 
can result in an uncertainty, due to the fact that HCN and HNC 
have slightly different dipole moments (3.0 vs. 3.3 D). 

Our density estimates are consistent with the estimates by 
Welch et al. ( 1 987), who estimate the cloud density to be >3 x 1 0^ 
cm"^ within the inner 1 pc and ~10'* cm"^ in the outer parts of 
the cloud. Our estimates are also consistent with (Plume et al. 
1997), who find - based on CS 2-1, 3-2, 5-4, 14-13, 10-9 tran- 
sitions - that «(H2) toward W49A varies between 5.2x10^ and 
1.8x10* cm"-'. Serabyn et al. (1993) estimate densities between 
2-6x10* cm"-^ for three clumps using CS (3-2, 5-4, 7-6, 10-9) 
and C-^'^S (5-4, 7-6, 10-9) transitions. Our mass, volume- and 
column density estimates also indicate a low filling factor for the 
dense gas. (Plume et al. 1997) probe masses and densities for 
a number of massive star-forming regions, and find, that filling 
factor for the dense gas is <25%. 

4.1. Possible heating sources 

The most probable mechanisms that contribute to the heating 
of the warm and dense molecular gas seen toward W49A are 
mechanical heating (shocks produced by stellar winds), and 
UV or X-ray irradiation. Another possible mechanism, that has 
been suggested for starburst galaxies, is heating by cosmic rays 
(Bradford et al. 2003). In the case of W49A, the only possible 
source of cosmic -rays above the Galactic background is the su- 
pernova remnant W49B, however, due to the high uncertainty in 
the distance of W49B (Brogan & Troland 2001), this remains 
an open question. Therefore, we focus on quantifying the effect 
of mechanical heating and irradiation by UV and X-rays for the 
warm and dense region revealed by our H2CO and HCN data. 

Mechanical heating The luminosity corresponding to the me- 
chanical heating produced by stellar winds of embedded O-type 
stars can be estimated from: 

L„ = 1.3xlO^*(— ^i^^ -)( ^-^ rfergs-i (3) 

llO-*MQyear-i/l2xl03kms-i/ ^ 

where Vw = 2000 km s"' can be applied as the velocity of the 
stellar wind (Tielens 2005). Peng et al. (2010) found evidence 
for two expanding shells with a size of ~2.9 pc, comparable to 
the size of the extended high excitation region we have found. 
They found that a constant mass-loss rate of ~ 1.2x10"* yr"' 
can sustain a wind-driven bubble with a size of ~2.9 pc, which 
corresponds to the mass loss rate of one O-type star. Using these 
values, the mechanical luminosity is = 1.56 x lO-'* erg s"' 
or about 405 Lq. Assuming an efficiency of 1-10 % (Loenen et 
al. 2008), the mechanical heating rate is ~4-40.5 Lq. However, 
this is only a lower limit of the luminosity corresponding to the 
mechanical heating. Alves & Homeier (2003) estimate that ~30 
O-type stars belong to the central stellar cluster, which has a 
diameter of ~6 pc. Taking this into account, the mass-loss rate, as 
well as the luminosity is about an order of magnitude larger than 
the value we derive by taking a mass-loss rate that belongs to 
one O-type star This gives 40.5-405 Lq as a mechanical heating 
luminosity. 

In addition to the effect of stellar winds, protostellar out- 
flows give a contribution to the mechanical heating. The largest 
outflow is connected to the most luminous water-maser emis- 
sion in the Milky Way (Gwinn et al. 1992). It has been char- 
acterized by Scoville et al. (1986) and recently by Smith et al. 
(2009). Adopting an outflow dynamical time-scale of 10"* yr, 
that corresponds to an average ~25 km s"' outflow speed (Smith 
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et al. 2009), an energy of 4.4x10'*^ erg (Scoville et al. 1986, 
scaled to a distance of 11.4 kpc), the mechanical luminosity is 
^mech = 1.4x10-'^ erg s"' or about ~400 L©. Assuming the same 
efficiency as for the mechanical luminosity produced by the stel- 
lar wind, only the most powerful outflow produces a mechanical 
luminosity that is 10 % of the mechanical luminosity produced 
by the stellar winds of an embedded cluster of ~30 O-type stars. 
This estimate gives a lower limit on the mechanical luminosity 
produced by outflows. Additional kinetic energy injection may 
have been resulted by past outflows of the current O-type stars. 
These flows may still be dissipating energy. The existence of 
these flows may be probed using ALMA. 

An additional source of the mechanical heating may exist 
in the form of shocks driven by the expansion of the (UC)Hii 
regions into the molecular cloud. The corresponding heating rate 
is proportional to the kinetic energy injection rate and is given 
by: 



2 df 2 



(4) 



where A is the surface area of the shocks, p is the density of the 
medium into which the shock is moving, and V is the velocity 
of the shock. The surface area of the shocks can be estimated 
based on De Pree et al. (1997), who measure the sizes of the 
(UC)Hn regions, based on high-resolution (O'.'S) 3.6 cm data. We 
use an average radius of 0.05 pc to calculate the surface area of 
the shocks. Based on «(H2) = 10* cm""*, a shock velocity of 
10 km/s results in a heating rate of 1.6 x lO-'^ for a cluster of 
30 O-type stars. For a shock velocity of 5 km/s, the mechanical 
heating rate of the expanding (UC)Hn regions corresponding to 
the cluster of 30 O-type stars is ~ 2 x lO-'* erg/s - in the same 
order of magnitude as the mechanical heating produced by the 
stellar winds of the 30 O-type stars. 



dence, that a population of small, hot dust grains with T ~350 K 
also exists, based on an excess in the near-infrared. The heat- 
ing rate of the gas per unit volume can be estimated using 
(Hollenbach & McKee 1979, 1989): 



= 1.2 xlO--,W^)'Vl^ 

I 1000 j I fl^n 



(5) 



x [1 - 0.8 exp(-75/ry„)](rd - rki„)(erg cm"^ s"') 

where the minimum grain size is set at a^m - 10 A. Using 
T~150 K and «=10* cm"-', assuming a volume filling factor of 
1-10% for the hot gas, a gas-to-dust ratio of 100, and that the 
fraction of the warm (T ~350 K) dust is up to 1% of the total 
dust, gas-grain heating for a 3x3 pc region toward the center of 
W49A is expected to be ~10^'*-10-'^ erg s"', equivalent to ~3-30 

Based on these estimates, mechanical heating and UV heat- 
ing seem to be the most probable heating mechanisms, while 
iiTadiation by X-rays and gas-grain collisions are less proba- 
ble. The contribution of radiation to the excitation of W49A has 
been investigated earlier by Roberts et al. (201 1), based on HCN, 
HNC and HCO^ transitions from the SLS. By comparing the ob- 
served line intensity ratios to PDR (Photon-dominated Region) 
and XDR (X-ray Dominated Region) models of Meijerink et al. 
(2007), they find, that HCN/HCO+ line ratios ai-e consistent with 
PDR models, while HCN/HNC line ratios are consistent with 
XDR models. They interpret this result as evidence that irradia- 
tion by UV and X-ray photons plays a minor role in the excita- 
tion of W49A. However, as seen above, based on the UV heating 
luminosity, UV radiation cannot be ruled out to contribute sig- 
nificantly to the heating of W49A. 



Radiative heating by the embedded stellar cluster The 

importance of heating by UV radiation from young massive 
stars, in the form of photo-electric emission by dust grains and 
PAHs, can be estimated from the far-infrared luminosity. Vastel 
et al. (2001) estimates a total FIR flux of 1.54 x lO"*" erg s"' 
cm"^, which, at the distance of W49A, is equivalent to a total IR 
luminosity of 6.3 x 10* Lq. Adopting an efficiency of 10"^ - 10"^ 
(Meijerink & Spaans 2005), the heating rate by UV irradiation 
is on the same order of magnitude as the luminosity that can be 
expected from the mechanical heating. 

X-ray heating Hard X-ray emission has been detected 
(Tsujimoto et al. 2006) toward the center of W49A, associ- 
ated with two Hii regions. They measure an X-ray luminosity 
of 3 X lO-'-' erg s"' in the 3-8 keV band. Assuming an efficiency 
of 10% (Meijerink & Spaans 2005), the X-ray heating luminos- 
ity is on the order of 0. 1 L©, about three orders of magnitudes 
below the luminosity of the mechanical and UV heating. In addi- 
tion, since X-ray emission is due to point-sources, rather than an 
extended region. X-rays mostly act locally and seem less likely 
to be responsible for the extended emission region what we de- 
tected. 

Gas-grain collisional heating At the densities estimated us- 
ing HCN transitions (10^ -10* cm""*), gas-grain collisions may 
give a significant contribution to the heating of the gas. Based on 
a greybody fit to the dust emission, Ward-Thompson & Robson 
(1990) derive a dust temperature of 50 K. They also find evi- 



4.2. Comparison to the cooling rate 

One of the main contributions to the total cooling in W49A is 
expected to occur via CO lines, therefore we derive a lower limit 
on the total cooling rate from the '-CO 3-2 (345.8 GHz) line 
luminosity based on data from the SLS. The CO line luminosity 
can be calculated from the CO 3-2 line intensities (Solomon & 
Vanden Bout 2005), as 

Leo [Kkms"' pc^] = 23.5Q[arcsec2](t^[Mpc])2/co[Kkms"'] (6) 
or alternatively, in units of L©: 

Leo [Lq] = 1.04X 10"^5eo[Jy]AV[kms-']v[GHz](4Mpc])^(7) 

where Ico = 593.4 K km/s, Tb = 34.5 K, Sco = 1.2 x 10'' Jy 
(using the Rayleigh-Jeans law to convert the brightness tem- 
perature to flux density) and AV = 17.2 km/s is the average 
line width. These numbers correspond to a spectnim smoothed 
to a FWHM = r beam, which covers the high-excitation re- 
gion around the center. Using these numbers, Leo(3-2) = 6.5 x 
10^ K kms ' pc^ = 9.65 Lq for the highly excited central square 
arcminute region. The fraction of the flux of the CO 3-2 tran- 
sition to the total CO line flux can be estimated using Radex, 
adopting an H2 volume density of 10* cm""* and a CO col- 
umn density of 10'^ cm"^. The CO column density estimate is 
based on the C'^O 3-2 transition (at ~337 GHz) and a ratio of 
'*O/'^O~1800 (Wilson & Rood 1994). For Lki,, = 100 K, the 
3-2 transition is expected to be 9%, and for Lkin = 300 K, 5 % 
of the total CO line flux. Using our Radex estimates of the total 
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CO line flux: 

Leo = 7.2 X lO'* - 1.3 X 10^ K kms"' pc' 
= 107.2 - 192.9 Lq 
= 4.1x10^^ -7.4x10^^ ergs"', 

gives a lower limit on the cooling, with an accuracy of ~20%, 
mainly as a contribution from the typical calibration error in de- 
riving the CO 3 -2 line intensity. 

Based on Neufeld et al. (2005), for kinetic temperatures 
>100 K and H2 volume densities ~10^ cm"^, the total cooling 
is up to an order of magnitude above the CO cooling, ~10^ Lq. 

Our estimates of the heating rate indicate that the two most 
important mechanisms to consider for W49A, are mechanical 
heating by the winds of O-type stars and UV-irradiation. These 
two mechanisms give a heating rate equivalent with 1000 Lq, 
depending on the heating efficiency. This is of a same order of 
magnitude as the expected total cooling, assuming that the total 
cooling rate is a factor of 10 higher than the CO cooling rate. 
This result supports that the effect of embedded O-type stars 
in the form of stellar winds and UV-irradiation can explain the 
high-excitation region around the center of W49A. 

4.3. Comparison to other warm and dense regions 

Galactic Center In our Galaxy, a similar level of excitation 
to what we found toward W49A was reported for clouds near 
the center. Huttemeister et al. (1993) measured kinetic temper- 
atures in 36 clouds near the Galactic center and based on NH3 
transitions, they have found a temperature component of ~200 
K, together with a cold gas and dust component. This sample, 
however, include no embedded stars, which could explain the 
highly excited gas component, and therefore, cannot be used as 
an analogue for W49A. De Vicente et al. (1997) have detected 
hot molecular gas toward the Sagittarius B2 molecular cloud. 
The highest temperatures measured toward Sagittarius B2 are in 
the range between 200-400 K, with typical densities of 10''- 10^ 
cm"^^, similar to what we find for the center of W49A. However, 
these conditions correspond to hot cores with sizes in the range 
of 0.5-0.7 pc, much smaller than the region with similar physi- 
cal properties in W49A. De Vicente et al. (1997) have also found 
a more extended component of warm (100-120 K) gas as a ring 
surrounding Sgr B2M and Sgr B2N with a radius of 2 pc and 
a thickness of 1.4 pc, with a density of about 2 x 10^ cm""*. 
Ferriere et al. (2007) study physical conditions in the innermost 
3 kpc of the Galaxy, and find a high-temperature gas component 
(T ~ 150 K), however, it corresponds to a low-density envi- 
ronment (nH2 ~ 10^'^ cm"^), unlike what we find for W49A. 
On Galactic scales, W49A shows a uniquely extended region of 
high (r > 100 K) temperature which appears to be tracing the 
feedback from the young and forming stars in the region on their 
surrounding gas. 

Starbursts and active galaxies In even nearby galaxies sen- 
sitivity and angular resolution mean that the physical condi- 
tions of the molecular gas can only be estimated on ~100 pc 
scales with uncertain filling factors. Nevertheless there are re- 
gions where temperatures similar to those in W49A are mea- 
sured. In particular, extended warm molecular gas (>100 K) has 
been detected in a number of starburst- and Seyfert galaxies (e.g. 
Rigopoulou et al. (2002)). Muhle et al. (2007) have detected 
a warm, Tkin ~ 200 K gas component with moderate density 
(n(H2) = 7 X lO-' cm"-') in the starburst galaxy M82, using p- 
H2CO line ratios, probing scales of ~320 pc. Ott et al. (2005) 



derived kinetic temperatures of 140 K and 200 K toward two po- 
sitions of the starburst galaxy NGC 253, using NH3(3,3) ATCA 
data on a scale of 65 pc. Mauersberger et al. (2003) have de- 
tected a warm (100-140 K) gas component toward the center 
of three nearby starburst galaxies (NGC 253, IC 342 and Maffei 
2) for regions of 300-660 pc, based on NH3 (l,l)-(6,6) tran- 
sitions. An even warmer component (T>400 K) was detected 
toward IC 342. Papadopoulos et al. (2011), based on a sample 
of 70 (Ultra-)Luminous Infrared Galaxies, most of them domi- 
nated by a warm (Tkin > 100 K) and dense (« > 10"* cm"-') gas 
phase with masses of ~10^ Mq, come to the conclusion that the 
main heating mechanism is rather cosmic-rays and highly super- 
sonic turbulence than far-UV radiation or shocks related to star- 
formation. This might be one of the key differences between the 
extended, but compared to scales of entire galaxies, still local 
phenomenon seen for W49A, and starburst galaxies. However, 
in the context of the Milky Way, W49A is a close analogue in 
terms of physical parameters to regions in starburst- and active 
galaxies. 

Another key comparison with warm (>100 K) and dense 
(>10^ cm"3) regions, such as starburst galaxies, is the initial 
mass function (IMF). Klessen et al. (2007) have found, based 
on hydrodynamical numerical calculations, that for regions of 
warm (~100 K) and dense (>10^ cm""*) gas, the slope of the 
IMF is in the range of -1.0 and -1.3, with a turn-over mass of 
7 Mq. Homeier & Alves (2005) derived an IMF in W49A with a 
slope of -1.70 + 0.30 and -1.6 + 0.3 down to masses of -10 Mq. 
This result may be consistent with the calculations of Klessen et 
al. (2007), and if so, supports that W49A is a Galactic template 
for starburst galaxies. However, incompleteness at low luminosi- 
ties and masses might explain the observed cut-off at masses of 

-10 Mq. 



Pressure as a diagnostic of starburst environments 

Environments with increased star-formation activity, such as 
starburst galaxies, massive star-forming regions and clouds near 
the Galactic center can be characterized with their pressures, 
P/k ~ iiT. Typical pressures that can be expected in the Galaxy 
and in the midplane of ordinary spirals are in the order of 10"* 
K cm"^ (Papadopoulos et al. 201 1). Table 4 shows the measured 
values of n, T and P/k for W49A, Sgr B2 and starburst galax- 
ies. P/k values show more than an order of magnitude variation. 
Based on our derived densities and temperatures, the gas pres- 
sure in W49A is up to 5.4x10** K cm"^, similar to those mea- 
sured toward SgrB2 and Arp 220. Other galactic nuclei appear 
to have lower pressures, which seems to be mainly the effect of 
a lower gas density. Beam dilution seems not to play a major ef- 
fect for the density estimates, given the high densities measured 
toward Arp 220 as well as the sensitivity, which makes measure- 
ments in external galaxies biased toward the warm and dense 
gas. However, the estimated densities significantly depend on 
the tracer used, as was shown by Greve et al. (2009), who find a 
10-100 factor difference between densities estimated from HCN 
compared to densities estimated from HCO^. Therefore, the un- 
certainty in P/k is about an order of magnitude. 
Pressures measured toward W49A and Sgr B2 are not entirely 
unique in our Galaxy. Dense gas has been detected toward a 
number of massive star forming regions, such as reported by 
Plume et al. (1997), Wu et al. (2010) and McCauley et al. 
(2011). However, the temperatures of those regions are only 
comparable to what we measure toward W49A locally', on sub- 
parsec scales, corresponding to hot cores, such as toward Orion 
KL (Mangum & Wootten 1993). Considering the spatial extent. 
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pressures measured toward W49A provide a Galactic analogue 
of starburst galaxies. 



Table 4. Comparison of physical parameters in W49A to 
Galactic center clouds and starburst galaxies 



Source 


n 


T 




P/k 






(K) 




(K 


W49A center 


1.8x10^ 


300 




5.4x10** 


W49A N clump 


5.6x10^ 


98 




5.5x10' 


W49A E tail 


1.2x10^ 


68 




8.2x10' 


W49A SW clump 


1.3x10' 


85 




1.1x10' 


Sgr B2" 


2x10' 


100 


-120 


2x10' -2.4x10' 


Henize 2-10* 


>lxlO'' 


50- 


100 


5 X 10' - 10^ 


M82'- 


10«_io« 


130 


-160 


10*"* - 10'' '^ 


Arp 220'' 


(O.S-DxlO" 


45- 


120 


1.4xl0'-1.2xl0'* 


Antennae"" 


10' 


90 




9x10^ 



" De Vicente et al. (1997) based on CHjCN J = 5-4, 8-7, 12-11 
lines. 

* Bayet et al. (2004), based on '-CO(J = 3-2, 4-3, 6-5, 7-6) and 
'^^CO(J = 3-2) 

Median values based on three different positions (Naylor et al. 
2010). Obtained from multiple species radiative transfer modeling of 
the lines of CS, HCO+, HCN, HNC, and C''^S. 

'' Greve et al. (2009), based on HCN and CS lines. 

' Zhu et al. (2003), based on CO and ''CO 2-1 and 3-2 lines. 



5. Summary and conclusions 

We have analyzed the physical conditions towards the central 
2x2 arcminutes region of W49A using JCMT maps with 15" 
resolution. We have reported a detection of warm (>100 K) and 
dense (>10^ cm"^) gas toward the center, extending up to a few 
parsecs. 

- We have characterized the excitation toward four main re- 
gions of W49A (center. Northern clump. South- West clump 
and Eastern tail) based on rotational diagrams of SO2, 
CH3OH and H2CO, taking account the optical depth and 
beam dilution. High excitation and column densities, as well 
as significant beam dilution is seen towards the center, likely 
to trace embedded hot cores. The off-center positions show 
excitation temperatures corresponding to about half of those 
derived for the center, and no significant beam dilution. 

- Based on 0-H2CO line ratios, we estimate kinetic temper- 
atures of the main regions of W49A. Toward the center, 
we derive a kinetic temperature of 250-300 K, 60-100 K 
toward the Northern clump, 80-150 K toward the South- 
west clump and 70-130 K toward the Eastern tail. Our es- 
timates show an approximately 1' x 1' (~3.3 x 3.3 pc) re- 
gion with extended warm (>100 K) gas towards the center of 
W49A. Based on HCN 3-2 / HCN 4-3 Hue ratios, the high- 
temperature central 1 ' x 1 ' region has H2 volume densities of 
>10^ cm"^*. The high excitation of W49A is comparable to 
clouds near the center of our Galaxy and to starburst galaxies 
with estimates on scales of ~100 pc. 

- The most impoitant heating mechanisms for W49A are me- 
chanical heating due to stellar winds of embedded, O-type 
stars and UV irradiation. X-ray irradiation and gas-grain 



collisional heating are less probable than other two mecha- 
nisms, while the effect of cosmic -rays remains an open ques- 
tion. However, the two main mechanisms corresponding to 
embedded young stellar populations provide a heating that 
no other mechanism is needed to contribute significantly in 
the heating of W49A. Based on earlier studies of starburst 
galaxies and their proposed heating mechanisms, W49A is a 
Galactic starburst analogue, in terms of excitation and phys- 
ical conditions, but not in terms of the heating mechanism 
which yields the high excitation, due to additional heating 
mechanisms that have been proposed to explain the excita- 
tion if starburst- and active galaxies (such as cosmic rays or 
supersonic turbulence). 

Our future work includes an analysis of the chemical in- 
ventory as well as kinematics of selected species, based on 
the entire SLS survey including data from the main survey 
(330-360 GHz, Plume et al. 2007) and will be presented in 
Nagy et al. (in prep). 
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Appendix A: Radex line ratio plots 

In sections 3.3 and 3.4 we compare the observed line ratios to 
line ratios calculated using the non-LTE code Radex (Van der 
Tak et al. 2007). Figure A. 1 shows the the 0-H2CO Sis -4i4 and 
533 - 432 line intensity ratios for H2 volume densities between 
10^ cm"-^ and 10^ cm"^ and temperatures between 20 K and 300 
K. Figure A. 2 shows the HCN 3-2 and 4-3 line intensity ra- 
tios for H2 column densities between lO-' cm""* and 10^ cm""* 
and temperatures between 20 K and 300 K. Both of these cal- 
culations use a molecular column density of 10'^ cm"^ and a 
background temperature of 2.73 K. We assume a line width of 
FWHM = 12 km s 'for H2CO , and a Une width of 18 km s 'for 
HCN for the central 40" x 40", in consistence with the observed 
line widths. 



H,CO 5„-4„ / H,CO 533-4 




Fig.A.l. Line ratios of 0-H2CO Sis - 4i4 and S33 - 432 as a 
function of kinetic temperature and H2 density for FWHM=12 
km s"'. The intensity scale is logarithmic. 
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Fig.A.l. Line ratios of HCN 3-2 and 4-3 as a function of ki- 
netic temperature and H2 density for FWHM=18 km s ' . 
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